BlenX-based compositional modeling of complex reaction mechanisms

Molecular interactions are wired in a fascinating way resulting in complex behavior of biological systems. Theoretical modeling provides a useful framework for understanding the dynamics and the function of such networks. The complexity of the biological networks calls for conceptual tools that manage the combinatorial explosion of the set of possible interactions. A suitable conceptual tool to attack complexity is compositionality, already successfully used in the process algebra field to model computer systems. We rely on the BlenX programming language, originated by the beta-binders process calculus, to specify and simulate high-level descriptions of biological circuits. The Gillespie's stochastic framework of BlenX requires the decomposition of phenomenological functions into basic elementary reactions. Systematic unpacking of complex reaction mechanisms into BlenX templates is shown in this study. The estimation/derivation of missing parameters and the challenges emerging from compositional model building in stochastic process algebras are discussed. A biological example on circadian clock is presented as a case study of BlenX compositionality

A genetic study based on PCNA-ubiquitin fusions reveals no requirement for PCNA polyubiquitylation in DNA damage tolerance

Post-translational modifications of Proliferating Cell Nuclear Antigen (PCNA) play a key role in regulating the bypass of DNA lesions during DNA replication. PCNA can be monoubiquitylated at lysine 164 by the RAD6-RAD18 ubiquitin ligase complex. Through this modification, PCNA can interact with low fidelity Y family DNA polymerases to promote translesion synthesis. Monoubiquitylated PCNA can be polyubiquitylated on lysine 63 of ubiquitin by a further ubiquitin-conjugating complex. This modification promotes a template switching bypass process in yeast, while its role in higher eukaryotes is less clear. We investigated the function of PCNA ubiquitylation using a PCNAK164R mutant DT40 chicken B lymphoblastoma cell line, which is hypersensitive to DNA damaging agents such as methyl methanesulfonate (MMS), cisplatin or ultraviolet radiation (UV) due to the loss of PCNA modifications. In the PCNAK164R mutant we also detected cell cycle arrest following UV treatment, a reduced rate of damage bypass through translesion DNA synthesis on synthetic UV photoproducts, and an increased rate of genomic mutagenesis following MMS treatment. PCNA-ubiquitin fusion proteins have been reported to mimic endogenous PCNA ubiquitylation. We found that the stable expression of a PCNAK164R-ubiquitin fusion protein fully or partially rescued the observed defects of the PCNAK164R mutant. The expression of a PCNAK164R-ubiquitinK63R fusion protein, on which the formation of lysine 63-linked polyubiquitin chains is not possible, similarly rescued the cell cycle arrest, DNA damage sensitivity, reduction of translesion synthesis and increase of MMS-induced genomic mutagenesis. Template switching bypass was not affected by the genetic elimination of PCNA polyubiquitylation, but it was reduced in the absence of the recombination proteins BRCA1 or XRCC3. Our study found no requirement for PCNA polyubiquitylation to protect cells from replication-stalling DNA damage. © 2017 Elsevier B.V

Cooperative interactions between p53 and NFkappaB enhance cell plasticity.

The p53 and NFkappaB sequence-specific transcription factors play crucial roles in cell proliferation and survival with critical, even if typically opposite, effects on cancer progression. To investigate a possible crosstalk between p53 and NFkappaB driven by chemotherapy-induced responses in the context of an inflammatory microenvironment, we performed a proof of concept study using MCF7 cells. Transcriptome analyses upon single or combined treatments with doxorubicin (Doxo, 1.5muM) and the NFkappaB inducer TNF-alpha (TNFalpha, 5ng/ml) revealed 432 up-regulated (log2 FC> 2), and 390 repressed genes (log2 FC< -2) for the Doxo+TNFalpha treatment. 239 up-regulated and 161 repressed genes were synergistically regulated by the double treatment. Annotation and pathway analyses of Doxo+TNFalpha selectively up-regulated genes indicated strong enrichment for cell migration terms. A panel of genes was examined by qPCR coupled to p53 activation by Doxo, 5-Fluoruracil and Nutlin-3a, or to p53 or NFkappaB inhibition. Transcriptome data were confirmed for 12 of 15 selected genes and seven (PLK3, LAMP3, ETV7, UNC5B, NTN1, DUSP5, SNAI1) were synergistically up-regulated after Doxo+TNFalpha and dependent both on p53 and NFkappaB. Migration assays consistently showed an increase in motility for MCF7 cells upon Doxo+TNFalpha. A signature of 29 Doxo+TNFalpha highly synergistic genes exhibited prognostic value for luminal breast cancer patients, with adverse outcome correlating with higher relative expression. We propose that the crosstalk between p53 and NFkappaB can lead to the activation of specific gene expression programs that may impact on cancer phenotypes and potentially modify the efficacy of cancer therapy

Loss of BRCA1 or BRCA2 markedly increases the rate of base substitution mutagenesis and has distinct effects on genomic deletions

Loss-of-function mutations in the BRCA1 and BRCA2 genes increase the risk of cancer. Owing to their function in homologous recombination repair, much research has focused on the unstable genomic phenotype of BRCA1/2 mutant cells manifest mainly as large-scale rearrangements. We used whole-genome sequencing of multiple isogenic chicken DT40 cell clones to precisely determine the consequences of BRCA1/2 loss on all types of genomic mutagenesis. Spontaneous base substitution mutation rates increased sevenfold upon the disruption of either BRCA1 or BRCA2, and the arising mutation spectra showed strong and specific correlation with a mutation signature associated with BRCA1/2 mutant tumours. To model endogenous alkylating damage, we determined the mutation spectrum caused by methyl methanesulfonate (MMS), and showed that MMS also induces more base substitution mutations in BRCA1/2-deficient cells. Spontaneously arising and MMS-induced insertion/deletion mutations and large rearrangements were also more common in BRCA1/2 mutant cells compared with the wild-type control. A difference in the short deletion phenotypes of BRCA1 and BRCA2 suggested distinct roles for the two proteins in the processing of DNA lesions, as BRCA2 mutants contained more short deletions, with a wider size distribution, which frequently showed microhomology near the breakpoints resembling repair by non-homologous end joining. An increased and prolonged gamma-H2AX signal in MMS-treated BRCA1/2 cells suggested an aberrant processing of stalled replication forks as the cause of increased mutagenesis. The high rate of base substitution mutagenesis demonstrated by our experiments is likely to significantly contribute to the oncogenic effect of the inactivation of BRCA1 or BRCA2.Oncogene advance online publication, 25 July 2016; doi:10.1038/onc.2016.243. © 2016 The Author(s

Minimum Criteria for DNA Damage-Induced Phase Advances in Circadian Rhythms

Robust oscillatory behaviors are common features of circadian and cell cycle rhythms. These cyclic processes, however, behave distinctively in terms of their periods and phases in response to external influences such as light, temperature, nutrients, etc. Nevertheless, several links have been found between these two oscillators. Cell division cycles gated by the circadian clock have been observed since the late 1950s. On the other hand, ionizing radiation (IR) treatments cause cells to undergo a DNA damage response, which leads to phase shifts (mostly advances) in circadian rhythms. Circadian gating of the cell cycle can be attributed to the cell cycle inhibitor kinase Wee1 (which is regulated by the heterodimeric circadian clock transcription factor, BMAL1/CLK), and possibly in conjunction with other cell cycle components that are known to be regulated by the circadian clock (i.e., c-Myc and cyclin D1). It has also been shown that DNA damage-induced activation of the cell cycle regulator, Chk2, leads to phosphorylation and destruction of a circadian clock component (i.e., PER1 in Mus or FRQ in Neurospora crassa). However, the molecular mechanism underlying how DNA damage causes predominantly phase advances in the circadian clock remains unknown. In order to address this question, we employ mathematical modeling to simulate different phase response curves (PRCs) from either dexamethasone (Dex) or IR treatment experiments. Dex is known to synchronize circadian rhythms in cell culture and may generate both phase advances and delays. We observe unique phase responses with minimum delays of the circadian clock upon DNA damage when two criteria are met: (1) existence of an autocatalytic positive feedback mechanism in addition to the time-delayed negative feedback loop in the clock system and (2) Chk2-dependent phosphorylation and degradation of PERs that are not bound to BMAL1/CLK

Compositional Modeling of Biological Systems

Molecular interactions are wired in a fascinating way resulting in complex behavior of bio-logical systems. Theoretical modeling provides us a useful framework for understanding the dynamics and the function of such networks. The complexity of the biological systems calls for conceptual tools that manage the combinatorial explosion of the set of possible interac-tions. A suitable conceptual tool to attack complexity is compositionality, already success-fully used in the process algebra field to model computer systems. We rely on the BlenX programming language, originated by the beta-binders process calculus, to specify and si-mulate high-level descriptions of biological circuits. Gillespie’s stochastic simulation algo-rithm applied for BlenX simulations requires the decomposition of phenomenological func-tions into basic elementary reactions. Systematic unpacking of complex reaction mechanisms into BlenX templates is shown. The estimation/derivation of missing parameters and the challenges emerging from compositional model building in stochastic process algebras are discussed. A biological example on circadian clock is presented as a case study of modeling

Neurospora crassa as a model organism to explore the interconnected network of the cell cycle and the circadian clock

Budding and fission yeast pioneered uncovering molecular mechanisms of eukaryotic cell division cycles. However, they do not possess canonical circadian clock machinery that regulates physiological processes with a period of about 24h. On the other hand, Neurospora crassa played a critical role in elucidating molecular mechanisms of circadian rhythms, but have not been utilized frequently for cell cycle studies. Recent findings demonstrate that there exists a conserved coupling between the cell cycle and the circadian clock from N.crassa to Mus musculus, which poses Neurospora as an ideal model organism to investigate molecular mechanisms and emerging behavior of the coupled network of the cell cycle and circadian rhythms. In this review, we briefly describe generic eukaryotic cell cycle regulation focusing on G1/S and G2/M transitions, and highlight that these transitions may be targeted for the circadian clock to influence timing of cell division cycle

Circadian rhythms synchronize mitosis in Neurospora crassa

The cell cycle and the circadian clock communicate with each other, resulting in circadian-gated cell division cycles. Alterations in this network may lead to diseases such as cancer. Therefore, it is critical to identify molecular components that connect these two oscillators. However, molecular mechanisms between the clock and the cell cycle remain largely unknown. A model filamentous fungus, Neurospora crassa, is a multinucleate system used to elucidate molecular mechanisms of circadian rhythms, but not used to investigate the molecular coupling between these two oscillators. In this report, we show that a conserved coupling between the circadian clock and the cell cycle exists via serine/threonine protein kinase-29 (STK-29), the Neurospora homolog of mammalian WEE1 kinase. Based on this finding, we established a mathematical model that predicts circadian oscillations of cell cycle components and circadian clock-dependent synchronized nuclear divisions. We experimentally demonstrate that G1 and G2 cyclins, CLN-1 and CLB-1, respectively, oscillate in a circadian manner with bioluminescence reporters. The oscillations of clb-1 and stk-29 gene expression are abolished in a circadian arrhythmic frq(ko) mutant. Additionally, we show the light-induced phase shifts of a core circadian component, frq, as well as the gene expression of the cell cycle components clb-1 and stk-29, which may alter the timing of divisions. We then used a histone hH1-GFP reporter to observe nuclear divisions over time, and show that a large number of nuclear divisions occur in the evening. Our findings demonstrate the circadian clock-dependent molecular dynamics of cell cycle components that result in synchronized nuclear divisions in Neurospora